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We present constraints on the nature of the first galaxies selected at 350 /xm. 
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350 /ill! (in the Bootes Deep Field) and also later serendipitous detections in 
the Lockman Hole. In determining multiwavelength identifications, the 350 fim 
position and map resolution of the second generation Submillimeter High An- 
gular Resolution Camera (SHARC II) are critical, especially in the cases where 
multiple radio sources exist and the 24 /im counterparts are unresolved. Spectral 
energy distribution templates are fit to identified counterparts, and the sample is 
found to comprise IR-luminous galaxies at 1 < 2; < 3 predominantly powered by 
star formation. The first spectrum of a 350 /im-selected galaxy provides an addi- 
tional confirmation, showing prominent dust grain features typically associated 
with star-forming galaxies. 

Compared to submillimeter galaxies selected at 850 and 1100 fim, galaxies se- 
lected at 350 /xm have a similar range of far-infrared color temperatures. However, 
no 350 /xm-selected sources are reliably detected at 850 or 1100 /xm. Galaxies in 
our sample with redshifts 1 < z < 2 show a tight correlation between the far- and 
mid-infrared flux densities, but galaxies at higher redshifts show a large disper- 
sion in their mid- to far-infrared colors. This implies a limit to which the mid-IR 
emission traces the far-IR emission in star-forming galaxies. 

The 350 yum flux densities (15 < 6*350 < 40mJy) place these objects near the 
Herschel/ SPIRE 350 /im confusion threshold, with the lower limit on the star 
formation rate density suggesting the bulk of the 350 fim contribution will come 
from less luminous infrared sources and normal galaxies. Therefore the nature 
of the dominant source of the 350 /im background — star-forming galaxies in the 
epoch of peak star formation in the universe — could be more effectively probed 
using ground-based instruments with their angular resolution and sensitivity of- 
fering significant advantages over space-based imaging. 

Subject headings: infrared: galaxies - submillimeter: galaxies - galaxies: star- 
burst - galaxies: high-redshift 



Introduction 



Submillimete r-selected galaxies (SMGs) were discovered in p i oneering lensed a nd blank- 

999h with the 



field surveys (e.g. , ISmail. Ivison &: Blainlll997l : iBarger et al.lll998l : lEales et al. 



850um-optimised Submillimeter Common User Bolometer Array (SCUBA; 



Holland et al. 



I999I] and later in similar surveys with mill imeter detectors (e.g., MAMBO, BOLOCAM, 



etc.: 



Bertoldi et al. 



galaxies (see, e.g.. 



2000 



Fox et al 



Laurent et al.ll2005l) . SM Gs mainly comprise ma ssive, star- forming 



2002) around 2 ~ 2 (IChapman et all l2005 l) with the bulk of 
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the emission generated at rest-frame far-IR wavelengths. SMGs are thus part of the IR- 
luminous galaxy population, which includes galaxies found in the local universe by the 



Infra r ed Astronornical S atellite {IRAS) All Sky Survey (ISoifer et al.lll984l : [Joseph fc Wright 



19851 : ISoifer et al.l 119871 ) and other galaxies detected in mid- and far-IR ba nds, most no 



tably (in terms of number selected) with the Infrared Space Observatory ( iKessler et al. 



19961). the Spitzer Spa ce Telescope (jWerner et a. 
jMurakami et al. 2007) out to ^ ~ 2 (see, e.g.. 



19991 : 
200i : 



Aussel. Elbaz. fc Cesarskv 



Le Floc'h et al. 



2004 



19991: 



Yan et al. 



2004r) and the Akari Infr a red Satellite 



Rowan- Robinson et al.l 119971: iPuget et al. 



Elbaz et al 



2004 



20021: 



Charv et al. 



Le Floc'h et al. 



200 



i 



200 



I 



Matsuhara 



lOnsdale et al. 



200 



The majority of SMGs have been selected at long submillimeter-millimeter bands (500- 
1300 /im). Shorter submillimeter wavelengths (200-500 /xm) are more demanding for ground- 
based observers. For example, on a good night at Mauna Kea, atrn ospheric transmission is 
about ~30% at 350 /xm but >80% at 850 /im; (ISerabyn et al.lll998l ). Despite this, the first 
galaxy selected purely by 350 /im emission — SMM J143206. 6 5-1-341613. 4 (=SSG 1, Short 
Submillimeter Galaxy 1) — was discovered (IKhan et al.l l2005l ) in a deep, blank survey of 
the Bootes Deep Fie ld with the second generation Submilli meter High Angular Resolution 
Camera (SHARC II; bowell et al. l boOsL iMoselev et al.ll2004r) along w ith a second detection, 
SMM J143206. 11+341648.4 (=SSG 2: Ixhanlbood : iKhan et aDl2007h . This survey, reaching 



13 mJy and currently the deepest at 350 /xm, obtained the first constraints on the 350 /xm 
source counts. Three additional 350 /xm-selected galaxies have been found in serendipi- 



(Laurent et al. 


2006; 


Coppin et al. 


2008) 



Short-wavelength submillim eter surveys are expected to principally select star-for ming 



galaxies at 1 < ^ < 3 (see, e.g., IPearson fc Khanll2009l: iKhan et al.l 120071: iKhanI |2006| ) . the 



epoch of peak star formation in the universe (e.g., iHopkins fc BeacomI 120061 ). Source count 
models that reproduce the observed 350 fim co unts include a predomi n antly IR-luminous 



galax y population evolving with redshift (e.g., iPearson fc KhanI l2009l : iFranceschini et al. 



20091 ) . but more observations are needed to verify these predictions. 



This paper presents a det ailed charac t erizat ion of each of the five 350 /im-selected galax- 
ies following the approach of iKhan et al.l (120051 ). Multiwavelength data are presented, and 
spectral energy distribution (SED) template fitting is used to provide constraints on the pho- 
tometric redshifts, thermal parameters (IR luminosity and dust temperature), and energy 
diagnostics. Additionally the first mid -IR spectruni of a 350 /xm-selected galaxy, obtained 
using Spitzer's Infrared Spectrograph (IHouck et al.l 120041 ). is given. The properties of the 
sample are compared with 850 and llOO/xm-selected SMGs. Table [1] identifies the five 
sources and includes short nicknames used for convenience. The WMAP first year cosmolog- 
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ical parameters (Hq = 71kms ^ Mpc ^, nni=0.27, Qx=0.73; Bennett et al. 2003) are used 
throughout this work. 



2. Observations, data reduction, and counterpart identification 

2.1. Observations 

The SHARC II observations, data reduction, and source extr a ction procedure for the 
Bootes Deep Field survey are discussed in de tail by iKhan et al.l (120071 ). Further details 
of source extraction are given by iKhanI (120061 ). Monte Carlo simulations for determining 
the survey completeness also provide a measure of the flux boosting. Approximately 4000 
artificial sources of random intensity and position were inse rted into the raw data map, then 
extracted using the same procedure as for the real sources (IKhan et al.ll2007l ). For recovered 
>3(T sources with input flux densities >15mJy, the ratio of measured flux density to input 
flux density indicates the average flux boosting is 1.07±0.23. Given the multiwavelength 
confirmation on SSG 1 and SSG 2, low outliers can be ruled out. Therefore we proceed 
without applying a correction factor to these sources and assume this also holds for the 
sources in the Lockman Hole. For the Lockman Hole, published 350 /xm flux densities were 
combined with public B, R, I, and z photometry from the Subaru Lock man Hole survey 



(iDye et al.l 120081 ). K band imaging from UKIDSS (iLawrence et a.1. 



20071) . archival Spitzer 



IRAC and MIPS imaging, and 1.4 GHz data from iBiggs &: IvisonI (120061 ). These compiled 
photometric data are presented in Table [H 

Additional 1.2 mm photometry of SSG 1 was obtained using MAMBO in 2005 January 
for ~3 hours in photometry mode. The horn antenna design produces incomplete sampling 
of the field, and because there was no jiggling, SSG 2 was not observed. The data were 
reduced using the standard package (mopsic). 

New VLA observations of the Bootes Deep Field were obtained on 2006 April 11. The 
VLA was in A configuration, giving a synthesized beam size of 1'.'36 x 1'.'50 at 1.5 GHz. Wide 
field imaging mode was used in order to avoid bandwidth smearing: two 25 MHz IFs of two 
polarizations each with 7 channels per IF. The phase center was located 10" north of SSG 1 
to avoid any possible phase center artifacts. Standard wide field imaging techniques were 
employed, including self-calibration and 3D corrections. The final image was generated with 
AIPS parameter robust = 1, giving an rms noise ~10 yuJy/beam. Automatic source search 
and flux density measurements were done with AIPS task sad. Both SSG 1 and 2 were 
detected, but neither was spatially resolved. 



Initial Spitzer/IRAC (IFazio et all 120041 ) and MIPS (IRieke et all l2004f ) data were ob- 
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tained during surveys of the NOAO Deep Wide Field in Bootes (lEisenhardt et al.ll2004j ). 
Significantly deeper observations were later obtained in 2006 February as part of the IRAC 
GTO program (Program ID 520). The IRAC observations consisted of six dithered 100 s 
frames in each IRAC field of view, covering the SHARC II survey area in all four IRAC 
bands. Multiplexer bleed and other detector artifacts were removed by applying the Spitzer 
Science Center's artifact mitigation code to the S14 version of the automated IRAC pipeline 
Basic Calibrated Data pr oducts. These cosmetically enhanc ed frames were then mosaiced 
using IRACProc v4.0beta (ISchuster. Marengo fc Pattenll2006l ) to ensure correct treatment of 
the noise for both point and extended sources. The MIPS observations consisted of eight 30 s 
cycles at 24 fim only and covered the entire 35 fim map area. The data were reduced follow- 
ing standard procedures (iGordon et al.ll2005l ). and source flux densities were measured with 
PSF fitting. The 99% confidence region of SSG 1 contains two visible-wavelength sources. 
Th e new 24 nm imagi ng shows that the source chosen as the most likely submm counterpart 
by iKhan et al.l (120051 ) is responsible for ~75% of the 24 /xm emission. 



The Spitzer/IRS observation of SSG 1 was designed according to the recommended SSG 
set-up. For each slit, there were six pointings along the slit at 24" spacing, all centered on 
the slit in its narrow dimension. This is equivalent to the normal point source method except 
that the target is observed in six different slit positions instead of the usual two, lessening the 
effects of fiat-fielding errors and bad pixels. In 3.1 hours (the total for the AOR, including 
peakups and overheads), there were five cycles of LLl (19.5-38.0 /im) and seven of LL2 
(14.0-21.3 yum) with 120 s ramps at each of the six slit positions. The data were reduced 



using the pipeline S15, cleaned using irsclean, and extracted with smart (IHigdon et al. 



20041 ). Two extraction techniques were compared: the first method involved extracting the 
spectra from the individual images and taking the median spectrum; the second method was 
to align the images, and then extract the spectrum. There was little difference in the final 
spectrum using either technique, and the spectrum is shown in Figure [TJ The most prominent 
features in the SSG 1 spectrum are due to polycyclic aromatic hydrocarbons (PAH). For 
quantitative analysis, PAH can be measured from integrating the flux above a baseline, but 
this method tends to underestimate the f lux density as corn pared to a profile decoraposit ion 
method such as the one used by pahfit ( ISmith et al.l 120071 : see also iGalliano et al.ll2008l for 
more detailed consideration of the two approaches). While both methods give similar PAH 
ratios and similar spatial variations of PAH strengths, absolute flux densities will differ due 
to assumptions about the PAH profiles. 
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2.2. Counterpart Identification 

The additional deep 24 fim and 1.4 GHz imaging of tlie Bootes field yi elded counterpart s 



in these bands for SSG 1 and 2 with low probabilities of chance association (iKhan et al.ll2007l ). 
The more precise positions were then used to identify optical and mid-IR counterparts. 
Finding counterparts for the three sources in the Lockman Hole is not straightforward: 
Figure [2] shows thumbnail images of the fields at various wavelengths. As can be seen, 
there exist multiple radio detections for two of the 350 fim objects. However the most likely 
counterparts can be determined with the following rationale: 



LH 350.1: the 1.4 GHz image shows two potential counterparts, but only the south- 
ern one is within the 9" SHARC H beam. In the visible images this appears to be 
a blend of two sources. The photometric redshift for this southerly source was ob- 
tained using HYPER-z (iBolzonella. Miralles fc PelFo I l2000f) in con i unctio n with the 
BRIzJK[3.6][4.5] photometry, similar to the approach of iDye et al.l (120081 ). The best- 
fitting solution is a starburst template at z = 1.19^q'H. The northern 1.4 GHz object 
is not associated with any visible-light detection, although it is prominent in the IRAC 
and MIPS images. Those images also show two more sources: an easterly object whose 
SED suggests it is a star or a low redshift elliptica l galaxy, and well to the west, the 
QSO RDS 054A aX z = 2.416 (jSchmidt et al.lll998l ). which is not expected to make a 
strong contribution to the 350 fim emission. 

LH 350.2: the 1.4 GHz image shows three sources within the beam. However, the 
two outermost objects are well-resolved in the visible imaging, suggesting that they 
are likely to be relatively nearby. Photo-z fitting puts both at 2; ~ 0.4, and SED 
template fitting puts a contribution to the 350 /im emission of <20% (assuming a dust 
temperature, > 30 K). Therefore, the very red middle object is likely to be the 
predominant source of the 350 /im emission with a best-fitting photo-z of 1.2lto2o- 

LH 350.3: there is a strong 24 /xm and weak 1.4 GHz source within 6" of the 350 /im 
position. The B through 4.5 fim images suggest that the mid-IR flux is a combination 
of two sources with the radio position lining up on the redder, western object. Best- 
fitting photometric redshifts are z = 2.47!^Q;g7 z = 0.56lo;3g for the western and 
eastern objects respectively. There remains the possibility of a merging system because 
the secondary solution for the eastern source is ^ = 2.60^o;33. Based on the radio 
identification, the western object is taken as the main source of the submillimeter 
emission. 
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3. Constraints on the photometric redshift and thermal emission 



The multiwavelength photometry in Table [2] was used to deri ve photometr i c red shifts 
and thermal emission parameters. Following the same approach as iKhan et al.l (120051 ) . the 
STARDUST2 template fit was used (jChania]||2009l ) to simultaneously obtain the photometric 
redshift ^phot, dust temperature T^, and the 8-1000 /im IR luminosity Ljr. Templates were 
based on local starbursts, and was minimised through a Levenberg-Marquardt technique. 
For SSG 1, the redshift obtained for the IRS spectrum in Figure [1] was used in the fit, 
but 2;spec = 1-05 agrees well with the 2;phot = 1-03 from the template fitting. The best-fit 
parameters are given in Table [31 and the best-fitting templates shown in Figure [21 

For SSG 1, SSG 2, and LH 350.3 the 350 fim flux density is well fit by the template, 
but in the cases of LH 350.2 and 3, the quality of the fit at 350 /im is clearly affected by 
the longer submillimeter data, which constrain the slope on the Ray leigh- Jeans side of the 
spectrum. Additionally, the 24 /xm flux density is useful in determining which side of the 
peak the 350 /xm datapoint lies. Usually the mid-IR is not used to constrain the restframe 
far-lR continuum emission because speciflc dust grains and AGN can have a major effect on 
mid-IR flux densiti es. Usage of the MIPS band was checked for a sample of SMGs (from 
Coppin et al.l l2008l ) by comparing template fits using measurements at 350 /im, 850 /im, 
and radio with measurements at 24 /im, 350 /im, and radio plus long-submillimeter limits. 
The two approaches agree within the stardust2 uncertainties. This indicate s our thermal 
parameters are not noticeably biased by the inclusion of the mid-IR data. (See IChania]||200g 
for more details.) 



For LH 350.3, the lack of a defining 1.6 /im rest-frame stellar continuum feature makes a 
significant AGN contribution plausible, although the relatively low radio fiux density would 
argue against this. Because of the absence of AGN templates in stardust2, the LH 350.3 
data were also fit us ing another routine with a mixture of AGN and starburst templates 
( jNegrello et al.ll2009l ). In this template fit, the favoured solution is a starburst at Zphot = 3.8 
(but without including the radio fiux density). Therefore we consider LH 350.3 to be, like the 
rest of this sample, predomina ntly star-forming. G iven the photometric redshift estimates 
from hyper-z, STARDUST2, and lNegrello et al.l (120091 ) . it is possible that LH 350.3 is the most 
distant galaxy so far discovered at short-submillimeter wavelengths. 



3.1. Energy source diagnostics 



As can be seen in Figure [31 all five sources are well fit in the mid-IR by templates that 
include features associated with star formation, i.e., a rest-frame 1.6 /im stellar continuum 
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bump and apparent PAH emission. For SSG 1, the presence of PAH is confirmed by the 
mid-IR spectrum, shown in Figure [H An emission hne from [Ne II] is hkely present though 
blended with possible PAH emission at 12.7 fim. The dust continuum is relatively flat until 
(rest wavelength) ~15.5/im. Although the 9.7/im silicate absorption is relatively weak, it 
is within the range seen in local IR-luminous galaxies. Table H] gives the emission feature 
fluxes, and a standard template (IRAS 22491 — 1808, an IR-luminous galaxy with moderately 
strong PAH emission) was used to obtain an overall best-fitting redshift from the spectral 
features of 1.05±0.01. 

The values of Lpah for the 7.7 and 11.3 fim r est-frame emission (6.8 and 2.0x10^ Lq) are 
consistent with the Lpah vs Ljr distribution of iPope et al. J2008h for a sample of 850 /im- 
selected SMGs. (This comparison used the integrated line flux densities with the Lir given 
in Table [3]) • The 7.7 /xm PAH luminosity also provide s an estimate of the star for mation rate 
(SFR) of 34OM0yr~^ following the relation given by lWeedman fc Houckl (120081 ). Given the 
strong star-forming features in the mid-IR spectrum, this comp ares favorably wit h the SFR 
derived from Ljr of 13QMqYt~^ following the relation given by iKennicuttI (119981 ). 



4. Discussion 



The different survey depths ma ke it complicated to produce fair statist ics for number 
counts. The Lockman Hole surveys (ILaurent et al.ll2006l : ICoppin et al.l 120081 ) were targeted 
at known sources, and sensitivity for serendipitous sources varied with distance from the 
target source. Even the blank fie ld survey in Bootes had sensitivity varying with location. 
(See Fig. 3 of iKhan et al.l 120071 .) However, simply considering the flux densities of the 
five detections and the total survey area (15.1 arcmin^)0, we estim ate resolving ~20% of 
the 350 fim contribution to the cosmic infrared background (GIB; iFixsen et al.l Il998l ) at 
5* > 17mJy. This fraction is reasonable given that the bulk of the background will be 
generated by galaxies near the break in the differential counts. It is also c onsistent with the 
30% estimate at 5* > 13mJy from the source counts of iKhan et al.l (120071 ). 



A lower limit on the star formation rate density (SFRD) can be found using the ther- 
mal parameters in Table [3] and the como ying volurne in 1 < z < 3. The conversion 
from far-IR luminosity to SFR is given by IKennicuttI (119981 ). The result for this epoch i s 
0.017Mq yr^^ Mpc~^, much lower than the values presented by iHopkins fc BeacomI (120061 ). 
This implies the bulk of the 350 /im contribution to the star formation rate density is from less 



5.8 arcmin for the Bootes observations; for LH 350.1 and 2 (jLaurent et al the map area is assu med 

to be that of the arguably similar SHADES observation of LH 350.3 of 3.1 arcmin^; ICoppin et al.ll2008l ) 
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luminous infrared galaxies and normal galaxies, consi stent with observations of IR-luminous 



galaxies in di fferent bands (e.g . . iLe Floc'h et al.ll2005[) and the predictions from source count 



models (e.g., Pearson &: Khan 20091 ) 



Our 350 /im-selected galaxies can be compared with sources selected at longer submil- 
limeter wavel e ngths . Suitable samples th at inc lude 350 and 2 4Atm photometry come from 
Coppin et al.l (120081 ) , iLaurent et al.l (120061 ) , and iKovacs et al.l (120061 ) . The photometric red- 
shifts and thermal parameters for these SMGs were again found using stardust2. Where 
there is no spect roscopic redshift, the photo-z was obt a ined from optical and mid-IR mag- 
nitudes gi ven bv iDve et al.l (120081) and lClements et al.l (120081 ). Following the Lj^-Td distri- 
bution by Blain et al. (j2004 ). there appear to be two trends separating the predominantly 
local IRAS galaxies from the higher redshift SMGs. (See Fig. |H) Compared with those 
selected at 850 and 1100 /im, the 350 /im-selected sources have similar dust temperatures 
and give no indication of temperature bias according to selection w avelength. Deeper 800- 
1000 /im imaging (e.g., with ALMA — IWootten fc Thompson! l2009h of the 350 /im-selected 
population would provide a more direct direct comparison of any selection biases. 

Mid-IR surveys to probe the 1 < ^ < 3 galaxy population have the advantage of bet- 
ter angular resolution, aiding counterpart identification (and hence more reliable redshift 
estimation), but a significant disadvantage is that without further far-IR observations, con- 
straints on the rest-frame far-IR thermal emission must be indirectly inferred. The 350:24 
flux density-redshift distribution for the SHARC II 350 /xm-selected SMGs compared with 
the 850 and 1100 /xm-selected SMGs is shown in Figure E] alongside colors from various tem- 
plates from lEfstathiou et al.l (120001 ). The overall distribution suggests that all samples are 
drawn from the same IR-luminous galaxy population, and the good agreement with the 
templates suggests these are predomtaantiy luminous and ultraluminous Infrared gala^esg 

Studies of infrared galaxies are often limited by unavailability of data at wavelengths 
longer than 24 /xm. In such cases, the IR luminosity can be es timated using only the 24 /^m 
flux density and various SED templates (jChary &: Elbazll200ll ) combined with an estimated 
redshift. For the SMG samples, the thermal emission estimates found this way compare 
surprisingly well with those from STARDUST2 for z < 1 as shown in Figure [Hi For these 
galaxies, the dominant errors on both luminosity predictors will arise from the photo-z 
fltting; these errors are highly non-Gaussian and non-trivial to determine. To some extent, 
the agreement in Figure [6] can be attributed to small number statistics and selection. The 
overall variation in the ratio over a wider redshift range could also imply a heterogenous SMG 
population. This would be consistent with the color in Figure but with a varying degree 



^10" < Lm < 10^2 L© and Lir > 10^^ L© respectively. 
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of AGN contribution. Another condition on this sample is the requirement of bright optical, 
mid-IR, and most importantly radio identifications; at high redshift, an AGN will boost the 
radio flux density and the emission in other bands (notably in the mid-IR), allowing this 
type of analysis to be performed. However none of the sources in the sample appears to be 
strongly affected by an AGN. At higher redshifts multiwavelength analysis becomes much 
more difficult. Some of the higher-^ 8 50/1100 /xm galaxies could simply have incorrectly 
identified 24 /im counterparts (see, e.g.. lYounger et al.ll2007l ). 



The restframe 12 /xm flu x density has been suggested to be a "pivot point" for SEDs 
(ISpinoglio et al.lll995l : see also lElbaz et al.ll2002l ). giving a fixed L12/I/FIR for both starbursts 
and AGN. This correlation was determined from local IRAS sources, but at z ~ 1, the 
rest frame 12 /im is shifted into the 24 fim band, and thus the relation at this redshift can 
be examined directly by our sample. We expect this relation to hold if the populations 
are similar. An observed difference in the populations is that distant mid-IR sources show 
stronger dust grain emission features around 12 /im compared to the lo cal IR-lurninous galaxy 
population. This difference could be due to differing AGN fractions (IWu et al.l 120091 ). 



Figure [7] shows the ratio of far-IR luminosity to observed Spitzer/MlPS 24 /im monochro- 
matic luminosity as a function of redshift. The ratio flattens out at 2; ~ 1 as the rest frame 
12 /im emission is shifted into the MIPS 24 fim band, but the dispersion increases at higher 
redshift. This suggests that the local empirical far-IR/12/im relation could apply to ear- 
lier epochs where the galaxy is undergoir ig rapid evolution, ra aking rest-frame 12 /im flux 
densities a better tracer of IR luminosity (ISpinoglio et al.lll995l ) than other mid-IR emission 
features, e.g., 7.7/im PAH selected at 2; ~ 2. At ^ ~ 1.4 the silicate absorption feature is 
redshifted into the MIPS 24 /im band and affects the flux density measurement. This effect 
is emphasized in Figure [5l where an excess in the 350/24 /xm colours is evident in the model 
SEDs, and galaxies in the redshift range where the silicate feature would cause a deficit in 
the 24 /im flux (1.3 < z < 1.6) are notably absent from Figure [71 

Various studies have sugges ted that the far-IR lu minosity can be deduced solely from 
mid-IR observations of galaxies. iBavouzet et al.l (120081 ) have shown that the far-IR luminos- 
ity of dusty galaxies observed with Spitzer correlate closely with the corresponding 8 and 
24 /im mid-IR luminosities. These studies were limited to redshifts z < 1 by the necessity 
to use Spitzer bands at 70 and 160 /im to calculate the far-IR luminosity. The mid-IR to 
far-IR luminosity relation was tentatively extended to higher redshifts {z ~ 2) by stacking 
the far-i nfrared fluxes in the long wavelength Spitzer bands, and these stacked results were 
used by ICaputi et al.l (120071 ) to calculate the bolometric luminosity function at 2; ~ 2. In 
contrast, our submillimeter observations directly access the high redshift Universe without 
the need for stacking images to obtain far-infrared fluxes. From Figure [7] we observe a tight 
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correlation between the far and mid-IR flux densities in our sample between 1 < z < 2, 
but there is a large dispersion in the infrared colours at z > 2. For z > 2, the 24 /xm band 
is sampling rest wavelengths shorter than 8 /im, and the far-IR luminosity cannot be accu- 
rately predicted from such observations. In order to confirm the mid-IR/far-IR relation for 
1 < 2; < 2 and to test the relation to higher redshifts, large samples of SMGs will be required. 
Such samples will be available following the large surveys planned with Herschel SPIRE (at 
250, 350, and 500 fxm: Icdfiin et al.ll20oi ) . 



5. Conclusion 

The discovery of the first 350 /zm-selected galaxies using SHARC II offers a different 
insight into IR-luminous galaxies than those selected in either the mid-IR or in longer sub- 
millimeter bands. The present sample comprises IR-luminous galaxies at 1 < ^ < 3 and 
resolves ~20% of the 350 /im background at S" > 17mJy. The implied lower limit on the star 
formation rate density suggests that the bulk of the 350 /im contribution is from galaxies of 
lower IR luminosity than the ones detected at present flux densities. 

The combination of 350 /im and mid-IR flux densities allows us to examine the relation 
of the mid and far-IR emission in dusty galaxies at high-redshift. There is surprisingly good 
agreement at z < 2, with the flattening of the IR:24/im/(l + z) luminosity ratio at z ^ 1 
evidence for the local empirical relation of rest frame 12 /im IR luminosity tracing IR lumi- 
nosity applying to distant IR-luminous galaxies. The wider dispersion at z > 2 suggests a 
limit to which the mid-IR can be used as a proxy for far-IR emission. We also find, by com- 
paring our 350 /im-selected sample to samples chosen at 850 /im and millimeter wavelengths, 
no evidence for a cold color temperature bias in SMGs. These results will be further tested 
through surveys with Herschel-SPIKE, which will detect SMGs in unprecedented numbers. 

The present 350 A^m-selecte d sample, with 15 < S^^^n < 40mJy, is near the S PIRE 350 /im 



blank- field detection limit (e.g.. |Pearson &: Kharul2009l : iFranceschini et al.ll2009l ) and also the 



detection limits of current longer-wavelength submillimeter instruments. Counterpart iden- 
tification is difficult even with the available angular resolution and will be far more difficult 
for space-based surveys with smaller telescopes. The achievable sensitivity and resolution 
of ground-based detectors will therefore provide a necessary complement to the space-based 
surveys, especially in characterising the nature of the faint SMG population. Whether the 
first 350 /im-selected galaxies are typical of the population will be probed through the combi- 
nation of surveys w ith SPIRE and future surveys from the ground (e.g., SCUBA 2-450 /im; 



Holland et al.ll2006f ). 
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Fig. 1. — The mid-IR spectrum of SSG 1 as observed by Spitzer IRS, indicated by black line. 
One-sigma error bars are shown. Wavelengths are in the rest frame based on the derived 
redshift z — 1.05. The spectrum of IRAS 22491—1808, used as a template to derive the 
redshift, is shown in grey. Prominent dust grain emission features are also indicated. The 
MIPS 24 iim passband is shown by the heavy black line below the spectrum. 
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Fig. 2. — Visible, mid-IR, and radio imaging of LH 350.1 (top), LH 350.2 (middle), and 
LH 350.3 (bottom). Each image is 20" on a side with a 9" diameter circle centered on the 
SHARC II 350 fira position. North is up and east to the left. Arrows in the LH 350.1 panels 
indicate the unrelated QSO RDS 054A [B and 20 cm panels), a star or low-redshift E galaxy 
{R panel), and the proposed SMG counterpart (4.5 /xm panel). Arrows in the LH 350.2 and 
350.3 panels indicate the proposed SMG counterparts. 
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Fig. 3. — Best-fitting SED templates from STARDUST2. Wavelengths are in the observed 
frame. 
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Fig. 4. — The IR luminosity-dust temperature distribution for various IR-luminous galaxy 
samples. For the 350 /im-selected galaxies, the red circles show SSG 1 and SSG 2, and 
triangles show LH 350.1, 350.2, and 350.3. Green stars represent the 850 and 1100 /im- 
selected SMGs for which th ere is 350 /xm fo l low-u p photometry. This is compared with 
850 /x m-selected SMGs f rom IChapman et al.l (120051] fopen blue cir c les) a nd IRAS galaxies 
from Dunne et al. ( 200ol ) (filled black squares) and Clements et al. ( 20091 ) (purple squares). 
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Fig. 5. — The 350/im:24/im flux density ratio as a function of redshift for various SMG 
samples as indicated in the legend. Lines show the same ratio for local SED templates from 
Efstathiou et~d1 tood ). 
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Fig. 6. — The ratio of best-fitting IR (8-1000 fim) luminosities obtained from the CE and 
STARDUST2 fitting, as a function of redshift. The diamonds denote the 350 /xm-selected 
galaxies, and filled circles denote the 850 and 1100 /xm-selected galaxies. The diamond with 
arrow denotes an upper hmit. Sources with redshifts beyond the useable range for the CE 
model (~3) are excluded from consideration. 
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Fig. 7. — Ratio of far-IR luminosity to monochromatic luminosity at observed wavelength 
24 /im as a function of redshift. Data are shown for the 350, 850, and 1100/im-selected 
SMGs. Numerator of the ratio is the rest frame far-IR (8-1000 fim) luminosity derived 
from the STARDUST2 template fitting. The denominator is the rest frame luminosity 
for the 24 /im observation, i.e., monochromatic luminosity uL^ at rest fr ame wavelength 
24 ij,m / {1 + z). At 2; ^ 1, observed 24 fim corresponds to rest 12 /im, which as lSpinoglio et al. 
( I1995I ) su ggested gives a g ood measure of total luminosity. The shaded region on the plot 
shows the ISpinoglio et al.l data for local luminous infrared galaxies as they would appear at 
z = 1. Filled triangles denote sources with spectroscopic redshift, while filled circles denote 
sources with photometric redshifts from the STARDUST2 template fit. 
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Tal)le 1: Tlie first five 350 //in-selectod galaxies 
Name [lAU] (nickname) Position [VLA] (J2000) 



Other names Position [350 /xm] ( J2000) 



SMM J143206.65+341613.4 (SSG 1) 


14 


32 


06.58 


+34 


16 


11.9 




14 


32 


06.65 


+34 


16 


13.4 


SMM J143206. 11+341648.4 (SSG 2) 


14 


32 


06.04 


+34 


16 


46.7 




14 


32 


06.11 


+34 


16 


48.4 


SMM J105308.3+571501 (LH 350.1) 


10 


53 


07.89 


+57 


15 


00.3 


SHARC II Source 3 


10 


53 


08.3 


+57 


15 


01 


SMM J105232.3+572448 (LH 350.2) 


10 


53 


32.26 


+57 


24 


47.4 


SHARC II Source 4 


10 


52 


32.3 


+57 


24 


48 


SMM J105243.2+572309 (LH 350.3) 


10 


52 


43.17 


+57 


23 


09.7 


LOCK 350.1 


10 


52 


43.2 


+57 


23 


09 



VLA positions are for identified VLA counterparts. The 350 //m position is used for the 

lAU name. 
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Table 2: Photometry for counterparts of 350 //m-selected galaxies 



' ^OOS 


SSG 1 


SSG 2 


LH 350.1 


LH 350.2 


LH 350.3 




<0.3 


1±0 02 


2 1±0 020 


07±0 007 


24±0 009 


R 


0.8±0.1 


0.3±0.1 


3 6±0 033 


27±0 013 


52±0 014 


I 


2.2 ± 0.1 


< 0.9 


5 5 + 051 


47 ± 017 


65 ± 018 


2, 






9 2 ± 085 


1 09 ± 040 


78 ± 036 


J 


10.5 ± 1.1 


< 7.9 


11.7 ± 0.65 


2.7 ± 0.45 


1.3 ± 0.41 


K 


42.0 ± 1.4 


< 11.5 


32.2 ± 0.89 


7.6 ± 0.56 


<1.5 


3.6 //m 


73.8 ± 0.3 


5.6 ± 0.3 


75.5 ± 7.0 


21.1 ± 2.0 


10.5 ± 1.1 


4.5 /im 


58.5±0.5 


7.4±0.5 


74.8 ± 6.9 


23.0 ± 2.2 


13.2 ± 1.5 


5.8 jim 


37.7±2.1 


8.1±2.1 


54.2 ± 7.0 


13.0 ± 4.5 


16.2 ± 5.2 


8.0 ^m. 


45.8±2.2 


3.0±2.2 


80.7 ± 8.6 


40.1 ± 6.0 


18.9 ± 5.6 


24 /im 


523.8±58.7 


187.2±29.4 


576 ± 13^^ 


204 ± 24^^ 


<80.2b 


70 yum [mJy] 


< 40 


< 40 








160 //m 


< 100 


< 100 








350 //m 


23.2±7.9 


17.1±6.4 


28.4±9.2 


37.0±13.4 


32.8±8.9 


850 //m 






<6 


<6 


<6 


1100//m 






1.1±1.4 


-0.2 ±1.4 


2.4±1.4 


1200 /im 


0.142±0.4 










20 cm [/iJy] 


49.8±11.4 


53.1±11.4 


52.5±5.2 


32.0±4.6 


15.8±4.8 



For non-detections, the flux density at the radio position is used. Upper limits are 3a. 
These data are used in the STARDUST2 SED template fitting. 



"Flux density is deblended using [3.6] /xm positions. 

''Flux density cannot be deblended, used as a limit instead. 
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Table 3: Best-fitting parameters from STARDUST2 template fitting. 



Nickname \og{LiR{SD2) / L^) [Lir{CE)Y Taust{SD2) Zphot{SD2) Ssm(SD2)'^ [mJy] 



SSG 1 

SSG 2 
LH 350.1 
LH 350.2 
LH 350.3 



11.9 



+0.1 
0.1 

i2.4i°:i 

1 9 q+0.10 

^-^•'-'-0.13 
-•-^•"-0.24 



[12.0] 

[12.6] 
[12.5] 
[11.9] 
<13.0] 



OQ 9+2.3 
OO. 2 3 

30.3l}-^ 
31.91^'^ 
34.11?:^ 



1 nQ+oi6c 

^■'-'•5-0.14 

1 07+0.8 

-0.3 



1.27^ 



-0.16 
-0.09 

1 26+°-^"' 

-■-•^"-0.12 
9 71 +0.03 



2.4±0.8 
2.7±0.7 
6.7±2.7 
2.4±0.6 
7.6±2.0 



"The luminosity from the Chary-Elbaz template fit, LmiCE), is obtained from the 24 /xm flux density and 
the STARDUST2 photometric redshift Zphot{SD2) (or spectroscopic redshift for SSG 1). 
''Predicted 850 /xm flux densities from the STARDUST2 fit. 
'^Spectroscopic redshift for SSG 1 is 1.05 ± 0.01. 
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Table 4: Features in the SSG 1 spectrum 



Feature Rest wavelength [/xm] Line Center [/i,m] Line Flux [10 ^^W cm 



PAH 

PAH 

[SIV] 

PAH 

PAH 

[Ne 11]^ 

[Ne V] 



7.7 

8.6 

10.5 

11.3 

12.7 

12.8 

14.3 



15.9 
17.4 
21.7 
23.1 

26.0 
26.9 
29.7 



14.5±4.4 
12.4±1.8 
1.2±1.8 
46.1±5.3 
3.0±0.6 
3.3±0.4 
0.5±0.3 



"The [Ne II] emission line could be blended with the 12.7 /um PAH feature. 



